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We study trapping and propagation of a matter-wave soliton through the interface between uni- 
form medium and a nonlinear optical lattice (NOL). Different regimes for transmission of a broad 
and a narrow soliton are investigated. Reflections and transmissions of solitons are predicted as 
function of the lattice phase. The existence of a threshold in the amplitude of the nonlinear opti- 
cal lattice, separating the transmission and reflection regimes, is verified. The localized nonlinear 
surface state, corresponding to the soliton trapped by the interface, is found. Variational approach 
predictions are confirmed by numerical simulations for the original Gross-Pitaevskii equation with 
nonlinear periodic potentials. 
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I. INTRODUCTION 



Investigation of processes of reflection, transmission 
and trapping of a nonlinear wave packet at the interface 
between two different nonlinear media represent one of 
the fundamental problems of the nonlinear physics [l|, . 
Recently, the problem of reflection and/or transmission 
of a soliton at the interface between a nonlinear uniform 
media and a linear periodic structure, under the con- 
ditions of the Bragg resonance, has been considered in 
Ref. There it was shown the possibility of controlling 
such a structure in the regime of the soliton mirror. 

Considering the actual experimental possibilities, it 
will be also of interest to study the transmission and 
trapping phenomena when we have a periodic variation 
in space of the parameter related to the nonlinearity. 
Such system can be realized in Bose-Einstein condensates 
(BEC) by using a periodically modulated in space ex- 
ternal magnetic field or optically induced Feshbach reso- 
nances Standing optical wave can induce in BEC 
periodic space modulation in the atom-atom scattering 
length. In the Gross-Pitaevskii (GP) equation it leads to 
periodic space modulations of the mean field nonlinear- 
ity; i.e., producing a NOL 7, 8, 9, 10, H, [H, [H, [11, [H . 
Considering the two-component case in the ID limit, the 
properties of BEC confined in NOL, as well as existence 
of soliton solutions and their stability, are investigated in 
Ref. 16]. BEC with finite segment of periodically space- 
modulated atomic scattering length (shallow optical lat- 
tice) are considered in Ref. [13 , where matter- wave opti- 
cal limits and bistability are predicted. For the dynamics 
of matter wave propagation under different conditions, 
see the review [ISj and references therein. Gap solitons 
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are analyzed in Ref. [l9| , where it was shown that local- 
ized nonlinear wave packets can exist in NOL for attrac- 
tive condensates (bright solitons) as well as for repulsive 
ones (dark solitons). The stability analysis showed that 
the bright solitons are stable in a very narrow region of 
parameters [20| . However, such analysis is absent in case 
of the existence of an interface. Linear surface states in 
lattices with management of the diffraction have been re- 
cently considered in Ref. [2l|. Surface soliton formation 
at an interface between two periodic media is studied in 
Ref. [H. 

In the present paper we consider regimes of reflec- 
tion, transmission and trapping of a matter wave soli- 
ton incident on the interface between uniform medium 
and a nonlinear optical lattice. Particular attention will 
be devoted to the possible existence of nonlinear surface 
states for matter waves. An interface induces changes in 
the effective potentials for the soliton center and width 
and can create a surface soliton. Also the stability can 
be enhanced. The dynamics of a BEC in a quasi-one- 
dimensional elongated trap will be treated by considering 
the GP formalism reduced to the one-dimensional (ID) 
space limit. 

Recently, an investigation done in Ref. [11] considered 
two-dimensional nonlinear surface states (surface soli- 
tons) at an interface in a superposition of a periodic po- 
tential and periodic modulations of the nonlinear space 
parameter. The physical system is motivated by optical 
structures writing on quartz by femtosecond laser (fs- 
laser). In this case, variations of the Kerr nonlinearity 
remain of the same sign and are out of the phase with 
the periodic variations of the linear refractive index. As 
opposed to this nonlinear optical system, in BEC case we 
can also realize the cases of periodic modulations when 
the nonlinearity changes sign. 

The transmission characteristics of solitons are de- 
fined by the effective potential induced by the interface 
and nonlinear periodic lattice. The effective potential 
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strongly depends on the soliton parameters. Unlike the 
soliton transmission through the linear lattice, in the case 
of nonlinear periodic lattice we have nontrivial intensity- 
dependence (number of atoms) for the form of the po- 
tential relief as well as a threshold behavior depending on 
the amplitude of the nonlinearity modulations in space. 
It means that by change in the amplitude of modula- 
tions performed by variation of external magnetic field 
near the Feshbach resonance point, we can form a mirror 
for the matter-wave solitons, selecting the solitons by the 
number of atoms. 

The paper is organized as follows. In section 2 the 
model is formulated and variational equations for the soli- 
ton parameters are derived. In section 3 characteristics of 
stationary soliton trapped by interface are investigated. 
The reflection and transmission regimes for narrow and 
broad solitons are analyzed. The summary of obtained 
results is given in conclusion. 



II. THE MODEL 

In order to describe the propagation of a matter wave 
soliton in the elongated quasi- ID condensate with attrac- 
tive interaction, we consider the GP equation in a ID 
space approach where the physical space-time variables 
are given by (x,t), and the corresponding dimensionless 
variables are {x,t): 

where ip = 'ip{x,t) and gioix) = 2fi,as(x)w^. Here, 
uj± is the transverse frequency of the trap and as{x) 
the spatially dependent atomic scattering length, which 
is supposed to vary in space for x > as as{x) — 
flo + 9{x) {So + aism{2kx)). Here, d{x) = 0(1) for x < 
0{x > 0), (5o is a constant, oq is the natural two-body 
scattering length, and the wave number k is related to 
the lattice period L by k = 2tt/L. The number of atoms 
N normalizes the wave-function as 



N = 



\iP\'^dx. 



(2) 



To avoid the collapse in the attractive BEG, the condition 
\as\N/uj± < 0.676 should be satisfied jH]. 

The transformation to the new set of dimensionless 
space-time variables (x, t) is given by the following: 

J — - / \ {x) ^ _ Sq 

X=kx, t = UJRt, j{x) = j-, Aq = -j r, 

|ao| |ao| 



Er = TvjJr,= — — , u = u{x,t) 
2m 



2uj± 



ao\^p. 



With the above, where Er is the recoil energy, we obtain 
the dimensionless form of the ID GP equation: 

iut + u^x + l{x)\u\'^u = 0, (3) 
j{x) = [70 + 9{x){Ao + 71 sin(2a;))], 



where 70 = ao/|ao| = ±1 (for the attractive and repulsive 
condensates respectively) and 71 = ai/|ao|. In the above 
equation and in the following, we use the abbreviated 
notation for partial differential equations, such that ut = 
du/dt. The normalization of u, N, relates to the number 
of atoms N, which is conserved. From Eqs. ^ and ([3]), 
we obtain: 



N = 



hk 



(4) 



Below we will consider the evolution of bright solitons 
(70 = 1). When solitons collide at the interface {x — 0), 
different scenarios are possible resulting in reflection, 
transmission or trapping. Let us consider different lim- 
iting cases of broad and narrow solitons (with respect to 
the period of modulations). To study the soliton evolu- 
tion we shall use the variational approach .25j . Accord- 
ing to this method we should calculate an averaged La- 
grangian and then, using the Euler-Lagrange equations, 
obtain the equations for the soliton parameters. 

The Lagrangian density corresponding to Eq. ([3]) is 
given by: 

L='-{utu*~utu)~K\' + ^j{x)\u\\ (5) 

In deriving our variational model we proceed from the 
following anzatz for a soliton: 



\/2ylsech 



(6) 



In order to obtain the equations for the soliton pa- 
rameters {A, a, (3, K, ^, (j)), we calculate the averaged La- 
grangian L = L(x, t)(ia; with the above trial function 

iV7o ^ AoN ^^llN^f, . 

where 

Fo{^, a) ^ 1+ tanh f ^ ) - ^ tanh^ f ^ 
3 \a J 3 \a 

/oc 
sech^(z)sin(2za-|-2^)dz, (8) 

with the integration variable z = {x — ^) / a. The Euler- 
Lagrange equations lead to the following: 



at = 4a/3, = 2k, 
N d 

«t = ——(^0^0+71^^2), 



8a 9^ 
Pt = -^f3^ + 



+ 



3N d /A0F0 + J1F2 



An^a da 



(9) 
(10) 

(11) 
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Eliminating the parameter (3 from the equations, we get 
the following evolution equation for the width a: 



Oitt 



16 



3iV d 
TT^ da 



^qFq + 71 Fa 
a 



This equation can be rewritten as: 



att 



da 



(12) 



(13) 



where 



N 



[47o + 3(AoFo+7iF2)]. (14) 



In a similar way for the soliton center we get: 



where 



N 

V^ia,^) = ^—{AoFo+-fiF2 
4a 



(15) 



(16) 



Typical profiles of effective potentials Va{a,£^) and 
V^(a,^) for the case of narrow solitons are shown in 
Fig. [1] The stationary point given by ao and is ob- 
tained in a self-consistent manner, using Eqs. (fT4)) and 
P^ . It should be noted that when the soliton norm 
decreases (increase in the width), the amplitude of the 
potential decreases. Accurate approximative analyt- 
ical expressions for the effective potentials Va{a,£,) and 
Vj(a,^), valid in a wide range of variables a and ^, can 
be obtained using the asymptotic representation of the 
integral ^ for F2 (?,"): 



F2{^,a)=g{a)sech\C/a) + 



2(1 + Q;2)7rQ!sin(2^) 
3 sinh(7ra) 



(l-Htanh(e/a)), 



(17) 



where g{a) = F2{0,a). This expression is obtained 
by sewing two different approximations for the integral, 
vaHd at the ranges: —00 < ^ ^ —A, ^ 1, and 
a ^ ^ < 00. The resulting interpolating formulae de- 
scribes well the integral F2 in all regions of ^, including 
the region < ^ < a. The dependence of the factor g{a) 
on the soliton width is shown in Fig. [51 



III. RESULTS 

All calculations below are performed for a solitary 
matter wave with background nonlinearity 70 = 1 and 
Ao = 0. In all the simulations of the soliton transmission 
and reflection, the incident wave packets are taken in the 
standard soliton form, with the initial parameters given 
by A = N/4, ao = l/A, (3 = 0, and = 0. In the starting 
position, the soliton is in a homogeneous medium. 
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FIG. 1: Profiles of the effective potentials K«(a,^o) [as a 
function of tfie width a, shown in (a)] and V^(ao,0 [as a 
function of ^, shown in (b)], where values of qo and ^0 cor- 
respond to the stationary point of the set of equations (|15|l . 
(|12p . Vertical line in (b) indicates the interface. Two cases 
are depicted. A*" = 1.5 with ao = 1.213, ^0 ~ 0.668 (dotted 
line) and iV = 2.5 with ao = 0.674, = 0.773 (sohd line). 
Other parameters are 70 = 1, 71 = 1.2, Ao — 0. 



A. Stationary soliton 

A stationary solution, within a semi-infinite lattice, is 
given by expression 



Uq 



(x) = V2Asech 



ao 



(18) 



where the stationary values, for the soliton position ^ and 
width a, are obtained in a self-consistent manner, from 
the set of equations 



dV^(a,0 



dVaia,0 



da 



= 0. 



(19) 
(20) 
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FIG. 2: Factor g{a) versus the soliton width a. 

Fig. [3] depicts stationary values of the soliton width ao 
and its position versus the norm N. 

As seen from Fig. |3{b), the position of the soliton cen- 
ter shifts from zero to some fixed point value as the num- 
ber of atoms increases. In order to check the stationary 
solution, which was obtained above, we solve directly the 
Gross-Pitaeskii equation Eq. with the starting wave 
packet (fT5)) . The results of the full PDE computation at 
t = 60 are given in Fig. 31 The parameters we consider 
are N — 2.5, 70 = 1 and 71 — 1.2. In this case, the 
stationary point obtained from the variational equations 
corresponds to ao = 0.6083 and ^0 = 0.7730. By starting 
with this variational point, the results of the full PDE 
simulation evolves to ao = 0.6293 and ^0 = 0.7778 at 
t — 60. Establishment of the stationary value of the soli- 
ton width in this case is shown in Fig.[5]for the evolution 
of the wave packet. 

We have also calculated the oscillation frequencies 
at the stationary point of the potentials Va{a,^) and 
V{(a, ^) versus the NOL strength 71 for the case = 2.5 
(ao = 1.6). The results of our calculations in the frame 
of the variational equations (|12p . and PDE simu- 
lations are depicted in Fig. [51 One can see satisfactory 
agreement between variational and PDE results. 



B. Reflection and transmission of narrow solitons 

Reflection and transmission of solitons at the interface 
between a homogeneous medium and the optical lattice 
can be described by Eq. considering the soliton as a 
particle travelling under the effective potential V^(a,^). 
Then, the condition of reflection or transmission is de- 
termined by the potential barrier height (see Fig. [1]). In 
Fig. [71 we present simulations for transmission and reflec- 
tion of the soliton at the interface, above and below the 
threshold value of the nonlinearity strength 71. The re- 
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FIG. 3: Dependences of the width ao (a) and position ^0 (b) 
of the stationary points on the norm N for the case 71 = 1.2. 



suits obtained with ODE calculations, based on Eq. (fT5|) . 
are shown with dotted lines. The solid lines correspond 
to PDE simulations, obtained from the GP equation ([3]). 
The soliton parameters are ao = 0.25, iV = 16 and v = 1, 
where u = = 2k; is the soliton velocity (see Eq. ([9])). 
The potential barrier heights, 71 = 0.022, indicated in 
Fig.jTjJa) and 71 — 0.04, indicated in Fig.jTfb), are respec- 
tively below and above the kinetic energy of the soliton, 
Ekin — v'^/2. As seen, the transmission and reflection 
conditions obtained in this case are in a good agreement 
with the PDE simulations of Eq. ([3]). We also observe 
that the true soliton dynamics deviates from the analyt- 
ical prediction for quite large times. At some depth of 
the soliton penetration into the optical lattice, the soliton 
can be trapped due to the radiation effects. 

To study numerically the travelling of a narrow soli- 
ton in a wide range of strengths of the optical lattice 
(0.022< 71 < 0.4) we carried out corresponding PDE 
simulations presented in Fig. ([8]). The soliton parame- 
ters were a = 0.25, A^ = 16, and v — 1. One can observe 
two regions of the optical lattice strength 71 which pro- 
vide trapping of the soliton: 71 < 0.023 and 71 > 0.335. 
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FIG. 4: The established stationary wave packet profile 114(2;) | 
a.t t = 60. Solid line stands for the full PDE results and 
squares are for the single soliton solution given by Eq. ((Gjl. 
The parameters are — 2.5, 71 — 1.2. The stationary point 
corresponds to ao = 0.6083 and ^0 ~ 0.7730 for the varia- 
tional approach; and ao = 0.6293 and ^0 ~ 0.7778 for the full 
PDE calculation at t = QO. 
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FIG. 5: Evolution of the wave packet width in the course 
of establishment of stationarity when starting from the sin- 
gle soliton solution given by Eq. The parameters are 
N = 2.5, ao = 0.6083, 71 = 1.2. 



The first region (71 < 0.023) corresponds to the soliton 
motion above the barrier. The trapping in this case is 
caused by an unavoidable radiation, which decreases the 
soliton kinetic energy, in the course of its motion in the 
optical lattice. The cause of the soliton trapping in the 
second case (71 > 0.335) can be explained by rearrange- 
ments of the wave packet due to deepening of the effective 
potential well, accompanied by strong radiation that re- 
sults in transformation of the incident moving soliton to 




FIG. 6: Frequencies cJa and uj^ versus the potential strength 
71. With solid lines we show results obtained with the varia- 
tional approach. The full PDE results are shown with squares. 
The soliton norm is A*' = 2.5. 

the stationary one. Typical profile of the trapped soliton 
in this case is depicted in Fig. The evolution of the 
soliton wave packet profile in the course of transmission 
of the optical lattice is presented in Fig. pU)) . As the 
soliton penetrates in the optical lattice, we observe that 
the amplitude of the transmitted soliton decreases due to 
noticeable radiation. 



C. Broad soliton 

Let us describe the dynamics of the broad soliton by 
Eqs. ([ini) and The VA is works well for the sohton 

propagation with the width less or of the order of the 
lattice period. The validity of VA for the dynamics of a 
broad (with respect to the lattice period) soliton should 
be checked by direct numerical simulations. As a rule 
one can expect a good agreement for the VA if radiation 
effects in propagation of the soliton in periodically mod- 
ulated media are small [2jj] . As shown in Ref . [9] , radia- 
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FIG. 7: Reflection and transmission of a soUton below (a) and 
above (b) the threshold value of 71 = 0.023. Solid lines stand 
for PDE simulations and dotted lines for variational ODE 
calculations. The parameters of the soliton are ao ~ 0.25, 

N = W,v = l. 
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FIG. 8: Travelling of the position of a narrow soliton at dif- 
ferent strengths of the periodic potential. Processes of trans- 
mission, reflection and trapping are depicted. The soliton 
parameters are ao = 0.25, A*' = 16 and v = 1. 




FIG. 9: Trapped soliton a,i t — 10. The simulation is started 
from the single soliton solution given by Eq. ([6} moving with 
velocity v = 1. Its initial position is ^0 = —5. Other parame- 
ters are A'' = 16, 70 = 1, 71 = 0.4, and ao = 0.25. 
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FIG. 10; Transmission of the narrow wave packet when 
starting from the single soliton solution with the parameters 
iV = 2.721, Qo = 1.47, ,^0 = -10, v = 1. The nonlinear 
optical lattice strength is 71 — 0.1. 



tivc effects at motion in the media with spatially periodic 
nonlinearity are small at the propagation of soliton with 
small velocity. The approximated expressions for the ef- 
fective potentials Va{a,£,) and V^{a,^) can be simplified 
when considering large a. So, as the soliton width grows 
to large a, the second term of Eq. ([T7]) can be neglected. 
By also imposing Aq = 0, the effective potential V^{a,£,) 
takes the form 

V^{a,0^^^^sech\^/a). (21) 

For 7i > 0, the potential V^{a, £,) is a potential well and, 
for 71 < 0, a potential barrier. It means that the reflec- 
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tion of the soliton becomes only possible provided that 
7i < 0. The sign of 71 is defined by the phase S of the 
periodic modulation of the nonlinearity ~ sin(2a; + 6). 
Thus, with the variation of such phase one can switch 
the matter-wave soliton from a transmission regime to a 
reflection one. The transmission (reflection) of the soliton 
occurs when the soliton kinetic energy is greater (smaller) 
than the potential barrier height. The threshold kinetic 
energy, Ecr, is given by 

Ecr^Y^ |Vc(a,e)|™a. = ^l7i|^- (22) 

Fig. pT|) depicts transmission of a broad soliton for the 
case 71 > 0, when the effective potential V(^{a, £,) is a po- 
tential well; and reflection and transmission of the soli- 
ton for the case 71 < 0, when the effective potential is 
a barrier. As seen, the conditions of transmission and 
reflection are well described by Eq. (|^^ . It should be 



The density is n « 10^m~^. The healing length and 
speed of sound are ^ « 2/im and c ~ 5mm/s. In a typical 
experiment [27j we could consider a soliton with about 
10'^ atoms and width « 2^ 4/im. In experiments we 
can vary the scattering length by using the Feshbach res- 
onance method by varying in space the external magnetic 
field B{x) near the resonant value Be, such that 

where ab is the background scattering length and A is 
the resonance width [28| . Other way to vary the scatter- 
ing length is the using of the optically induced Feshbach 
resonances [1, @ . Typical values of variations of the 
external field B are in the interval (0 — 200)/im with the 
period 10/im around the value 352 G, where the scatter- 
ing length has the minimal value w — 0.23nm. 
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FIG. 11: Numerical simulations of reflection and transmis- 
sion of a broad soliton with the parameters A'^ — 0.6 and 
Qfo = 6.667. The solid line stands for the case 71 > 0, when 
the efi^ective potential V^{a, ^) is a potential well. Dotted and 
dashed lines are for the case 71 < 0, when the effective poten- 
tial is a barrier for an incident soliton. The soliton velocities 
i; = 2k = 0.075 (dotted line) and v = 2k = 0.05 (dashed line) 
correspond to the soliton kinetic energies upper and below the 
potential barrier, respectively. 

noted that trapping of the soliton is not observed since 
the effective potential V^{a,^) is a short-range one, such 
that, far from the interface, motion of the soliton can be 
considered as free. 

Let us discuss the possible effects which can be pre- 
dicted in experiments with BEG in a cigar-type trap. 
As an example, we can consider the ''Li condensate in 
the elongated trap with the transverse frequency uj± « 
27r X lO'^Hz and the longitudinal frequency uJx ~ few Hz. 



IV. CONCLUSION 

We have investigated reflection, transmission and trap- 
ping of a matter-wave soliton, which is propagating 
through the interface between a nonlinear uniform me- 
dia and a nonlinear optical lattice. We study analytically 
two different limits for broad and narrow soliton dynam- 
ics. In both the cases we obtain the characteristics of a 
soliton trapped by an interface, corresponding to a lo- 
calized nonlinear surface wave. We derive the effective 
potentials for the soliton center-of-mass and the width 
induced by the joint action of the interface and the me- 
dia periodic nonlinearity. We obtain the parameters of a 
localized nonlinear surface state, corresponding to a soli- 
ton trapped near the interface. Near the stationary point 
of effective potentials Eq. (fH)) and Eq. (|16p . we have cal- 
culated frequencies of oscillations of the trapped soliton 
center and its width. It was also obtained the thresh- 
old value of the NOL strength 71, separating the trans- 
mission and reflection regimes for incident solitons. The 
predicted surface soliton states can be observed in exper- 
iments with BEG in optically induced NOL and in non- 
linear optical systems with periodically modulated Kerr 
nonlinearity. 
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